Proximity of the Superconducting Dome and the Quantum Critical Point in the 
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We use the dynamical cluster approximation to understand the proximity of the superconducting 
dome to the quantum critical point in the two-dimensional Hubbard model. In a BCS formalism, 
T c may be enhanced through an increase in the d-wave pairing interaction (Vd) or the bare pairing 
susceptibility (xod)- At optimal doping, where Vd is revealed to be featureless, we find a power-law 
behavior of Xod(w = 0), replacing the BCS log, and strongly enhanced T c . We suggest experiments 
to verify our predictions. 



Introduction- The unusually high superconducting 
transition temperature of the cuprates remains an un- 
solved puzzle, despite more than two decades of intense 
theoretical and experimental research. Central to the ef- 
forts to unravel this mystery is the idea that the high 
critical temperature is due to the presence of a quan- 
tum critical point (QCP) which is hidden under the su- 
perconducting dome Pi. Numerical calculations in the 
Hubbard model, which is accepted as the de-facto model 
for the cuprates, strongly support the case of a finite- 
doping QCP separating the low-doping region, found to 
be a non- Fermi liquid (NFL) , from a higher doping Fermi- 
liquid (FL) region [2, 3]. Calculations also show that in 
the vicinity of the QCP, and for a wide range of tem- 
peratures, the doping and temperature dependence of 
the single-particle properties, such as the quasi-particle 
weight [2j, as well as thermodynamic properties such as 
the chemical potential and the entropy, are consistent 
with marginal Fermi liquid (MFL) behavior jij. This 
QCP emerges by tuning the temperature of a second- 
order critical point of charge separation transitions to 
zero and is therefore intimately connected to q = charge 
fluctuations [B| . Finally, the critical doping seems to be in 
close proximity to the optimal doping for superconduc- 
tivity as found both in the context of the Hubbard [5| 
and the t-J model Even though this proximity may 
serve as an indication that the QCP enhances pairing, 
the detailed mechanism is largely unknown. 

In this Letter, we attempt to differentiate between two 
incompatible scenarios for the role of the QCP in super- 
conductivity. The first scenario is the quantum critical 
BCS (QCBCS) formalism introduced by She and Zaanen 
(She-Zaanen) [7|. According to this, the presence of the 



QCP results in replacing the logarithmic divergence of 
the BCS pairing bubble by an algebraic divergence. This 
leads to a stronger pairing instability and higher critical 
temperature compared to the BCS for the same pairing 
interactions. The second scenario suggests that remnant 
fluctuations around the QCP mediate the pairing inter- 
action [1, Q. In this case the strength of the pairing 
interaction would be strongly enhanced in the vicinity 
of the QCP, leading to the superconducting instability. 
Here, we find that near the QCP, the pairing interac- 
tion depends monotonically on the doping, but the bare 
pairing susceptibility acquires an algebraic dependence 
on the temperature, consistent with the first scenario. 

Formalism- In a conventional BCS superconductor, 
the superconducting transition temperature, T c , is deter- 
mined by the condition Vxq(w = 0) = 1, where x'o ^ s 
the real part of the q = bare pairing susceptibility, and 
V is the strength of the pairing interaction. The transi- 
tion is driven by the divergence of Xo( w = 0) which may 
be related to the imaginary part of the susceptibility via 
Xq(oj = 0) = - J dwxo(u)/w- And Xofa) itself can be 
related to the spectral function, A^{uj), through 
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where the summation of £ G { — 1, 
to anti-symmetrize Xo( w )- 111 a FL, 
AT(o;/2)tanh(a;/4T), and Xo( T ) « N(0)ln(uj D /T) with 
iV(0) the single-particle density of states at the Fermi 
surface and u>n the phonon Debye cutoff frequency. 
This yields the well known BCS equation T c — 
cj D exp[-l/(N(0)V)]. In the QCBCS formulation, the 
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BCS equation is Vx'(w = 0) = 1, where x' is fully dressed 
by both the self energy and vertices associated with the 
interaction responsible for the QCP, but not by the pair- 
ing interaction V. In the Hubbard model the Coulomb 
interaction is responsible for both the QCP and the pair- 
ing, so this deconstruction is not possible. Thus, we will 
use the more common BCS T c condition to analyze our 
results with Vxb(w = 0) = 1 where xb 1S dressed by 
the self energy but without vertex corrections. Since the 
QCP is associated with MFL behavior, we do not expect 
the bare bubble to display a FL logarithm divergence. 
Here, we explore the possibility that Xo( w = 0) ~ l/T a . 
The two-dimensional Hubbard model is expressed as: 
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where cl (cfc ff ) is the creation (annihilation) operator for 
electrons of wavevector k and spin a, n; L<y — cJ CT Cj CT is 
the number operator, e° = — 2t (cos(k x ) + cos(k y )) with 
t being the hopping amplitude between nearest-neighbor 
sites, and U is the on-site Coulomb repulsion. 

We employ the dynamical cluster approximation 
(DCA) [10( to study this model with a Quantum Monte 
Carlo (QMC) algorithm as the cluster solver. The DCA 
is a cluster mean-field theory which maps the original lat- 
tice onto a periodic cluster of size N c = L 2 embedded in 
a self-consistent host. Spatial correlations up to a range 
L c are treated explicitly, while those at longer length 
scales are described at the mean-field level. However the 
correlations in time, essential for quantum criticality, are 
treated explicitly for all cluster sizes. To solve the cluster 
problem we use the Hirsch-Fye QMC method [HI HH and 
employ the maximum entropy method (l3j to calculate 
the real-frequency spectra. 

We evaluate the results starting from the Bethe- 
Salpeter equation in the pairing channel: 
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where x is the dynamical susceptibility, Xo{Q)p [= 
—G(P+Q)G(—P)] is the bare susceptibility, which is con- 
structed from G, the dressed one-particle Green's func- 
tion, r is the vertex function, and indices pt—l and ex- 
ternal index Q denote both momentum and frequency. 
The instability of the Bethe-Salpeter equation is detected 
by solving the eigenvalue equation Txo4> = A<^> [H| for 
fixed Q. By decreasing the temperature, the leading A 
increases to one at a temperature T c where the system 
undergoes a phase transition. To identify which part, xo 
or r, dominates at the phase transition, we project them 
onto the d-wave pairing channel (which was found to be 
dominant For Xo, we apply the d-wave projection 

as XodM = EfeXo(w,g = 0) k g d (k) 2 / 9d(k) 2 , where 
gd(k) = (cos(fc^) — cos(fcj,)) is the d-wave form factor. 
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FIG. 1: (Color online) Plots of leading eigenvalues for differ- 
ent channels at the critical doping for N c — 12 and N c — 16 
site clusters. The inset shows the phase diagram with su- 
perconducting dome, pseudogap T* and FL Tx temperatures 
from Ref. @] 



As for the pairing strength, we employ the projection as 
V d = J2k,k' 9d(k)T k ^g d (k')/y ik g d (k) 2 , using T at the 
lowest Mastsubara frequency [l6| . 

To further explore the different contributions to the 
pairing vertex, we employ the formally exact parquet 
uations to decompose it into different components [161 
Namely, the fully irreducible vertex A, the charge 
(S=0) particle-hole contribution, $ c , and the spin (S=l) 
particle-hole contribution, through: T = A + $ c + $ s . 
Similar to the previous expression, one can write V d — 
+ V d + VT* , where each term is the d-wave component 
of the corresponding term. Using this scheme, we will be 
able to identify which component contributes the most 
to the d-wave pairing interaction. 

Results- We use the BCS-like approximation, dis- 
cussed above, to study the proximity of the supercon- 
ducting dome to the QCP. We take U = 6t (At = 1) on 
12 and 16 site clusters large enough to see strong evidence 
for a QCP near doping 5 ps 0.15 0, 13, Q- We explore the 
physics down to T s» 0.1 1J on the 16 site cluster and 
T « 0.07J on the 12-site cluster, where J 0.11 [H| is 
the antiferromagnetic exchange energy. The fermion sign 
problem prevents access to lower T. 

Fig. [TJ displays the eigenvalues of different channels 
(pair, charge, magnetic) at the QC filling. The results for 
the two cluster sizes are nearly identical, and the pairing 
channel eigenvalue approaches one at low T, indicating a 
superconducting d-wave transition at roughly T c = 0.007. 
However, in contrast to what was found previously (l6j |. 
the q = charge eigenvalue is also strongly enhanced, 
particularly for the larger N c = 16 cluster, as it is ex- 
pected from a QCP emerging as the terminus of a line of 
second-order critical points of charge separation transi- 
tions Q. The inset shows the phase diagram, including 
the superconducting dome and the pseudogap T* and FL 
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FIG. 2: (Color online) Plots of Vd, the strength of the d-wave 
pairing interaction for various temperatures with U = 1.5 
(At = 1) and N c = 16. Vd decreases monotonically with 
doping, and shows no feature at the critical doping. In the 
inset are plots of the contributions to Vd from the charge 
V^ and spin V£ cross channels and from the fully irreducible 
vertex Vd versus T at the critical doping. As the temperature 
is lowered, T <C J ~ 0.11, the contribution to the pairing 
interaction from the spin channel is clearly dominant. 



Tx temperatures. 

In Fig. [21 we show the strength of the d-wave pair- 
ing vertex Vd versus doping for a range of temperatures. 
Consistent with previous studies [l9[ , we find that Vd falls 
monotonically with increasing doping. At the critical 
doping, 6 C = 0.15, Vd shows no feature, invalidating the 
second scenario described above. The different compo- 
nents of Vd at the critical doping versus temperature are 
shown in the inset of Fig. [2J As the QCP is approached, 
the pairing originates predominantly from the spin chan- 
nel. This is similar to the result of Ref. [16[ where the 
pairing interaction was studied away from quantum crit- 
icality. 

In contrast, the bare d-wave pairing susceptibility xod 
exhibits significantly different features near and away 
from the QCP. As shown in Fig. [3l in the underdoped 
region (typically 8 = 0.05), the bare d-wave pairing sus- 
ceptibility x'odi 1 ^ = 0) saturates at low temperatures. 
However, at the critical doping, it diverges quickly with 
decreasing temperature, roughly following the power-law 
behavior 1/vT, while in the overdoped or FL region it 
displays a log divergence. 

To better understand the temperature-dependence 
of Xodi 10 — 0) a t t ne QC doping, we looked into 
T 1 - 5 Xo' d (w)/w and plotted it versus uj/T in Fig.@] When 
scaled this way, the curves from different temperatures 
fall on each other such that T x - 5 ^ d {u)/u) = H(uj/T) m 
(w/T)- 15 for uj/T >9« At/ J. For < uj/T < At/ J, the 
curves deviate from the scaling function H(x) and show 
nearly BCS behavior, with Xod(w)/w| Ci;=0 which is weakly 
sublinear in 1/T as shown in the inset. The curves away 
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FIG. 3: (Color online) Plots of x'odi^ — 0), the real part of 
the bare d-wave pairing susceptibility, at zero frequency vs. 
temperature at three characteristic dopings. The solid lines 
are fits to x'odi^ = 0) = B/VT + A ln(u c /T) for T < J. In 
the underdoped case (5 = 0.05), x'odi 1 ^ = 0) does not grow 
with decreasing temperature. At the critical doping (5 = 8 C = 
0.15), Xodi^ = 0) shows power-law behavior with B = 0.04 
for the 12 site, and B — 0.09 for the 16-site clusters (in both 
A — 1.04 and ui c — 0.5). In the overdoped region (8 — 0.25), 
a log divergence is found, with B — obtained from the fit. 
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FIG. 4: (Color online) Plots of r 15 Xo' d (w)/a; versus lj/T at 
the QC doping (8 = 0.15) for JV C = 16. The arrow denotes the 
direction of decreasing temperature. The curves coincide for 
uj/T > 9 ~ (At/ J) defining a scaling function H(ui/T), cor- 
responding to a contribution to Xod(T) = ~ J cL)Xod( w )/ w 
1/VT as found in Fig. [3] For u/T > 9 w (At/ J), H(u/T) w 
(lo/T)~ 1,5 (dashed line). On the x-axis, we add the label 
T 3 /T w (At/ J), where T 3 represents the energy scale where 
curves start deviating from H. The inset shows the un- 
sealed zero-frequency result Xo'<i( w )/ a; L = o plotted versus in- 
verse temperature. 



from the critical doping (not displayed) do not show such 
a collapse. In the underdoped region (8 = 0.05) at low 
frequencies, x'odi^ )/ 00 S oes to zero with decreasing tem- 
perature (inset). In the FL region (8 = 0.25) Xod{ ljJ )/ ul 
develops a narrow peak at low ui of width uj ss Tx and 
height tx 1/T as shown in the inset. 
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Discussion- x'odi^)/ 10 reveals details about how the 
instability takes place. The overlapping curves found 
at the QC filling contribute a term T-^Hfo/T) to 
Xo d (w)/w or x' d( T ) K 1 /^T as found in Fig. O There 
is also a component which does not scale, especially at 
low frequencies. In fact, Xod{ u} )l u a ^ zer0 frequency in- 
creases more slowly than 1/T as expected for a FL. From 
this sublinear character, we infer that the contribution 
of the non-scaling part of x 'd( w )/ w to the divergence of 
Xo d (T) is weaker than BCS and may cause us to overes- 
timate A and underestimate B in the fits performed at 
the critical doping in Fig. [3j In addition, if H(0) is finite, 
it would contribute a term to Xod(-^) that, increases like 
1/T 1 - 5 , so H{0) = 0. From Eq. [T] we see that the con- 
tribution to XQ d (u;)/aj at small co comes only from states 
near the Fermi surface. H(0) = would indicate that 
the enhanced pairing associated with x'od(T) °^ 1/VT is 
due to higher energy states. The vanishing of Xod( w )/ a; 
in the pseudogap region (8 — 0.05) for small frequency 
when T — > indicates that around the Fermi surface, the 
dressed particles do not respond to a pair field. Or, per- 
haps more correctly, none are available for pairing due 
to the pseudogap depletion of electron states around the 
Fermi surface. Thus, even the strong <i-wave interaction, 
seen in Fig. [3J is unable to drive the system into a su- 
perconducting phase. In the overdoped region, x'odi^/^ 
displays conventional FL behavior for T < Tx, and the 
vanishing Vd suppresses T c . 

Together, the results for xod and Vd shed light on the 
shape of the superconducting dome in the phase diagram 
found previously [f|. With increasing doping, the pair- 
ing vertex Vd falls monotonically. On the other hand, 
Xod(T) is strongly suppressed in the low doping or pseu- 
dogap region and enhanced at the critical and higher dop- 
ing. These facts alone could lead to a superconducting 
dome. Futhermore, the additional algebraic divergence 
of xbd(T) seen in Fig. [3] causes the superconductivity to 
be enhanced even more strongly near the QCP where one 
might expect T c cx (V^B) 2 , with B = — J dxH(x), com- 
pared to the conventional BCS form in the FL region. 

Similar to the scenario for cuprate superconductivity 
suggested by Castellani et al. 0] , we find that the super- 
conducting dome is due to charge fluctuations adjacent 
to the QCP related to charge ordering. However, we dif- 
fer in that we find the pairing in this region is due to an 
algebraic temperature dependence of the bare suscepti- 
bility xod rather than an enhanced ci-wave pairing vertex 
V d , and that this pairing interaction is dominated by the 
spin channel. 

Our observation in the Hubbard model offers an exper- 
imental accessible variant of She-Zaanen's QCBCS. We 
use the bare pairing susceptibility xo while She-Zaanen 
use the full x> which includes all the effects of quan- 
tum criticality but not the correction from the pairing 
vertex (the pairing glue is added separately). This de- 
composition is not possible in numerical calculations or 



experiments since both quantum criticality and pairing 
originate from the Coulomb interaction. However, the 
effect of quantum criticality already shows up in the one- 
particle quantities, and the spectra have different be- 
haviors for the three regions around the superconduct- 
ing dome. She-Zaanen assume that x"( w ) 1/^ Q for 
T s < u) < lo c , where ui c is an upper cutoff, and that 
it is irrelevant (a < 0), marginal (a — 0), or relevant 
(a > 0), respectively in the pseudo gap region, FL re- 
gion and QCP vincity. We find the same behavior in xo 
and we have the further observation that near the QCP 
T s « (At/J)T and a = 0.5. 

Experiments combining angle-resolved photo emission 
(ARPES) and inverse photo emission results, with an 
energy resolution of roughly J, could be used to construct 
Xod and explore power law scaling at the critical doping. 
Since the energy resolution of ARPES is much better 
than inverse photo emission, it is also interesting to study 
Xg d (w)/w| w= o! which only requires ARPES data, but not 
inverse photo emission. 

Conclusion- Using the DCA, we investigate the d- 
wave pairing instability in the two-dimensional Hubbard 
model near critical doping. We find that the pairing in- 
teraction remains dominated by the spin channel and is 
not enhanced near the critical doping. However, we find a 
power-law divergence of the bare pairing susceptibility at 
the critical doping, replacing the conventional BCS log- 
arithmic behavior. We interpret this behavior by study- 
ing the dynamic bare pairing susceptibility which has a 
part that scales like x ' d (u;)/u; ~ T~ 1 - b H(uj /T), where 
H(lu/T) is a universal function. Apparently, the NFL 
character of the QCP yields an electronic system that is 
far more susceptible to d-wave pairing than the FL and 
pseudogap regions. We also suggest possible experimen- 
tal approaches to exploit this interesting behavior. 
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